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Abstract 

In the Standard Model, the coupling of the Higgs boson to b quarks is weak, 
leading to small cross sections for producing a Higgs boson in association with 
b quarks. However, Higgs bosons with enhanced couplings to b quarks, such 
as occur in supersymmetric models for large values of tan [5, will be copiously 
produced at both the Tevatron and the LHC in association with b quarks which 
will be an important discovery channel. We investigate the connections be- 
tween the production channels, bg — > bh and gg — > bbh, at next-to-leading 
order (NLO) in perturbative QCD and present results for the case with two 
high-pT b jets and with one high-pr 6 jet at both the Tevatron and the LHC. 
Finally, the total cross sections without cuts are compared between gg — > bbh 
at NLO and bb -> h at NNLO. 

1. Introduction 

In the Standard Model, the production of a Higgs boson in association with b quarks is suppressed by 
the small size of the Yukawa coupling, g^bh = m^/v ~ 0.02. However, in a supersymmetric theory with 
a large value of tan j3, the 6-quark Yukawa coupling can be strongly enhanced, and Higgs production in 
association with b quarks becomes the dominant production mechanism. 

In a four-flavor-number scheme with no b quarks in the initial state, the lowest order processes are 
the tree level contributions gg — > bbh and qq — > bbh, illustrated in Fig. 1. The inclusive cross section for 
gg — > bbh develops potentially large logarithms proportional to L5 = log(Q 2 / m b) which arise from the 
splitting of gluons into bb pairs. 1 Since Q ^> nib, the splitting is intrinsically of 0{a s L\ ) ), and because 
the logarithm is potentially large, the convergence of the perturbative expansion may be poor. The 
convergence can be improved by summing the collinear logarithms to all orders in perturbation theory 
through the use of b quark parton distributions (the five-flavor-number scheme) [4] at the factorization 
scale fip = Q. This approach is based on the approximation that the outgoing b quarks are at small 
transverse momentum. Thus the incoming b partons are given zero transverse momentum at leading 
order, and acquire transverse momentum at higher order. In the five-flavor-number scheme, the counting 
of perturbation theory involves both a s and l/Lj. In this scheme, the lowest order inclusive process is 
bb — > h, see Fig. 2. The first order corrections contain the 0(a s ) corrections to bb — > h and the tree 

'it should be noted that the b mass in the argument of the logarithm arises from collinear bb configurations, while the large 
scale Q stems from b transverse momenta of this order, up to which factorization is valid. The scale Q is the end of the collinear 
region, which is expected to be of the order of [1, 2, 3]. 



level process gb — > bh, see Fig. 3, which is suppressed by 0(1/ Lb) relative to bb — ► h [5]. It is the 
latter process which imparts transverse momentum to the b quarks. The relevant production mechanism 
depends on the final state being observed. For inclusive Higgs production it is bb — ► h, while if one 
demands that at least one b quark be observed at high-py, the leading partonic process is gb — > 6/i. 
Finally, if two high-pj- 6 quarks are required, the leading subprocess is gg — > 

The leading order (LO) predictions for these processes have large uncertainties due to the strong 
dependence on the renormalization/factorization scales and also due to the scheme dependence of the b- 
quark mass in the Higgs 6-quark Yukawa coupling. The scale and scheme dependences are significantly 
reduced when higher-order QCD corrections are included. 

Section 2 describes the setup for our analysis, and in Section 3 we compare the LO and NLO QCD 
results for the production of a Higgs boson with two high-p^ b jets. Section 4 contains a discussion of the 
production of a Higgs boson plus one high-pr b jet at NLO, including a comparison of results within the 
four-flavor-number and the five-flavor-number schemes. We consider the corresponding inclusive Higgs 
cross sections in Section 5. Although motivated by the MSSM and the possibility for enhanced b quark 
Higgs boson couplings, all results presented here are for the Standard Model. To a very good approxima- 
tion the corresponding MSSM results can be obtained by rescaling the bottom Yukawa coupling [6, 7]. 
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Fig. 1: Sample Feynman diagrams for gg — > bbh and qq — > bbh production. 




Fig. 2: Feynman diagram for bb — > h production. 





Fig. 3: Feynman diagrams for gb — > bh production. 



2. Setup 

All results are obtained using the CTEQ6L1 parton distribution functions (PDFs) [8] for lowest order 
cross sections and CTEQ6M PDFs for NLO results. The top quark is decoupled from the running of 
rrib(fJ.) and a s (/j,) and the NLO (LO) cross sections are evaluated using the 2 (l)-loop evolution of a s (fi) 
with a^ LO (Mz) = 0.118. We use the MS running b quark mass, m,b(n), evaluated at 2 (l)-loop 
for (Jtvlo (°loX with the b pole mass taken as m& = 4.62 GeV. The dependence of the rates on the 
renormalization ([Ir) and factorization (hf) scales is investigated [5, 6, 7, 9, 10] in order to estimate the 
uncertainty of the predictions for the inclusive Higgs production channel and for the Higgs plus 1 6-jet 
channel. The dependence of the Higgs plus 2 b- jet rates on the renormalization ((J,r) and factorization 



(^f) scales has been investigated elsewhere [6, 7] and here we fix \x = /j,r = hf = (Zmb + Mh)/4, 
motivated by the studies in Refs. [1, 2, 3, 5, 6, 7, 9, 10]. 

In order to reproduce the experimental cuts as closely as possible for the case of Higgs plus 1 or 2 
high-pr b quarks, we require the final state b and b to have a pseudorapidity | r] \ < 2 for the Tevatron and 
| 77 |< 2.5 for the LHC. To better simulate the detector response, the gluon and the b/b quarks are treated 
as distinct particles only if the separation in the azimuthal angle-pseudorapidity plane is AR > 0.4. For 
smaller values of AR, the four-momentum vectors of the two particles are combined into an effective b/b 
quark momentum four- vector. All results presented in the four-flavor-number scheme have been obtained 
independently by two groups with good agreement [6, 7, 11, 12]. 



3. Higgs + 2 b Jet Production 

Requiring two high-py bottom quarks in the final state reduces the signal cross section with respect to 
that of the zero and one b-tag cases, but it also greatly reduces the background. It also ensures that 
the detected Higgs boson has been radiated off a b or b quark and the corresponding cross section is 
therefore unambiguously proportional to the square of the 6-quark Yukawa coupling at leading order, 
while at next-to-leading order this property is mildly violated by closed top-quark loops [6, 7]. The 
parton level processes relevant at lowest order are gg — > bbh and qq — > bbh, as illustrated in Fig. 1. 
Searches for the neutral MSSM Higgs bosons h, H, A produced in association with b quarks have been 
performed at the Tevatron [13]. 

The rate for Higgs plus 2 high-pr b jets has been computed at NLO QCD in Refs. [6, 7] and is 
shown in Fig. 4 for both the Tevatron and the LHC. The NLO QCD corrections modify the LO predictions 

by < 30% at the Tevatron and < 50% at the LHC. The total cross section plots include a cut on p^ b > 20 
GeV, which has a significant effect on the cross sections. We show the dependence of the cross section 
on this cut in Fig. 5. The NLO corrections are negative at large values of the cut on p^ b and tend to be 
positive at small values of p^ b . 
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Fig. 4: Total cross sections for pp(pp) — ► bbh + X at the Tevatron and the LHC as a function of the Higgs mass Mh with two 
high-pT & jets identified in the final state. The b/b quarks are required to satisfy p b J b > 20 GeV. We fix fj, — fj,R — (if — 
(2m b + M h )/4. 



4. Higgs + 1 b Jet Production 

The associated production of a Higgs boson plus a single b quark (or b quark) is a promising channel for 
Higgs production in models with enhanced bbh couplings. The cross section is an order of magnitude 
larger than that for Higgs plus 2 high-pr b jet production for the cuts imposed in our analysis. 
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Fig. 5: Total cross sections for pp(pp) — > bbh + X at the Tevatron and the LHC as a function of the cut pTcut in p^/ b for a 
Higgs mass M h — 120 GeV with two high-pr b jets identified in the final state. We fix fi = fiR = hf = (2ra(, + M h )/4. 



In the four-flavor-number scheme, this process has been computed to NLO, with the momentum 
of one of the b quarks integrated over [6, 11, 12]. This integration yields a potentially large factor 

L5. Both the total cross sections and the dependence on the p^ b cut at the Tevatron and the LHC are 
shown in Figs. 6 and 7. The NLO corrections increase the cross section by < 50% at the Tevatron 
and < 80% at the LHC. The renormalization/factorization scales are varied around the central value 
A* = = = (2mf, + M/j)/4. At the Tevatron, the upper bands of the curves for the four-flavor- 
number scheme in Figs. 6 and 7 correspond to [Ir = = 2/x, while the lower bands correspond to 
Vr = A*F = m/2. The scale dependence is more interesting at the LHC, where the upper bands are 
obtained with = n/2 and \ip = 2/j,, while the lower bands correspond to = 2/x and pp = jjl/2. 
At both the Tevatron and the LHC, the width of the error band below the central value (/j, = /j,r = /xp) 
is larger than above. 

In the five-flavor-number scheme, the NLO result consists of the lowest order process, bg — > bh, 
along with the 0(a s ) and 0(1/ Lf,) corrections, which are of moderate size for our scale choices [9]. 
The potentially large logarithms L5 arising in the four-flavor-number scheme have been summed to all 
orders in perturbation theory by the use of b quark PDFs. In the five-flavor-number scheme, the upper 
bands of the curves for the Tevatron in Figs. 6 and 7 correspond to /j,r = \i and fip = 2\i, while the lower 
bands correspond to //r = /u/2 and up = [i. At the LHC, the upper bands are obtained with \ir = \i 
and [ip = 2/x, while the lower bands correspond to [ir = 2/j, and fip = n/2. The two approaches agree 
within their scale uncertainties, but the five-flavor-number scheme tends to yield larger cross sections as 
can be inferred from Figs. 6 and 7. 

Contributions involving closed top-quark loops have not been included in the five-flavor-number 
scheme calculation of Ref. [9]. This contribution is negligible in the MSSM for large tan/3. In the 
four-flavor scheme, the closed top-quark loops have been included and in the Standard Model reduce the 
total cross section for the production of a Higgs boson plus a single b jet by ~ —7% at the Tevatron and 
13 % at the LHC for M h = 120 GeV [11, 12]. 

5. Inclusive Higgs Boson Production 

If the outgoing b quarks are not observed, then the dominant process for Higgs production in the five- 
flavor-number scheme at large values of tan (3 is bb — > h. This final state contains two spectator b quarks 
(from the gluon splittings) which tend to be at low transverse momentum. At the LHC this state can 
be identified through the decays into fi + fi~ and t + t _ for the heavy Higgs bosons H, A at large values 
of tan (5 in the MSSM [14]. The bb -> h process has been computed to NLO [5] and NNLO [10] in 
perturbative QCD. The rate depends on the choice of renormalization/factorization scale Hr/f, and at 
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Fig. 6: Total cross sections for pp(pp) — > bbh + X at the Tevatron and the LHC as a function of the Higgs mass Mh 
with one high-pT b jet identified in the final state. The b(b) quark is required to satisfy p b r/, b > 20 GeV. We vary the 
renormalization/factorization scales around the central value /j, = (j,r = /if = {2mt + Mh)/4 as described in the text. 
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Fig. 7: Total cross sections for pp(pp) — > 6Wi + X at the Tevatron and the LHC as a function of the cut pTcut in b for 
a Higgs mass Mh = 120 GeV with one high-pT b jet identified in the final state. We vary the renormalization/factorization 
scales around the central value fj, = [in = = (2rrt6 + Mh)/4 as described in the text. 



NLO a significant scale dependence remains. The scale dependence becomes insignificant at NNLO. It 
has been argued that the appropriate factorization scale choice is /j,p = (M^ + 2mb)/4 [2, 3] and it is 
interesting to note that at this scale, the NLO and NNLO results nearly coincide [10]. 

An alternative calculation is based on the processes gg — > bbh and qq — > bbh (four-flavor-number 
scheme), which has been calculated at NLO [6, 11, 12]. Despite the presence of the logarithms Lh 
in the calculation based on gg — > bbh, which are not resummed, it yields a reliable inclusive cross 
section, as evidenced by Fig. 8. A sizeable uncertainty due to the renormalization and factorization scale 
dependence remains which might reflect that the logarithms are not resummed in this approach, so 
that the perturbative convergence is worse than in the corresponding case of tih production [15]. In the 
Standard Model, the closed top-quark loops have been included in the four-flavor-number calculation and 
reduce the inclusive NLO total cross section for pp(pp) — > bbh by ~ —4% at the Tevatron and ~ —9% 
at the LHC for M h = 120 GeV [11, 12]. In the MSSM, the closed top quark loops are negligible for 
large tan/3 [6, 7]. 

The NLO four-flavor-number scheme calculation is compared with the NNLO calculation of 
bb — > h (five-flavor-number scheme) in Fig. 8. The two calculations agree within their respective scale 
uncertainties for small Higgs masses, while for large Higgs masses the five-flavor-number scheme tends 



to yield larger cross sections. Note that closed top-quark loops have not been included in the NNLO 
calculation of bb — > h [10]. 

To all orders in perturbation theory the four- and five-flavor number schemes are identical, but 
the way of ordering the perturbative expansion is different and the results do not match exactly at finite 
order. The quality of the approximations in the two calculational schemes is difficult to quantify, and 
the residual uncertainty of the predictions may not be fully reflected by the scale variation displayed in 
Fig. 8. 
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Fig. 8: Total cross sections for pp(pp) — » bbh + X at the Tevatron and the LHC as a function of the Higgs mass Mh with 
no 6 jet identified in the final state. The error bands correspond to varying the scale from /j,r = = (2mf, + Mh)/8 to 
fJ.R = (J.F = {2m b + M h )/2. The NNLO curves are from Ref. [10]. 



6. Conclusions 

We investigated bbh production at the Tevatron and the LHC, which is an important discovery channel 
for Higgs bosons at large values of tan j3 in the MSSM, where the bottom Yukawa coupling is strongly 
enhanced [13, 14]. Results for the cross sections with two tagged b jets have been presented at NLO 
including transverse-momentum and pseudorapidity cuts on the b jets which are close to the experimen- 
tal requirements. The NLO corrections modify the predictions by up to 50% and reduce the theoretical 
uncertainties significantly. For the cases of one and no tagged b jet in the final state we compared the 
results in the four- and five-flavor-number schemes. Due to the smallness of the b quark mass, large 
logarithms might arise from phase space integration in the four-flavor-number scheme, which are 
resummed in the five-flavor-number scheme by the introduction of evolved b partem densities. The five- 
flavor-number scheme is based on the approximation that the outgoing b quarks are at small transverse 
momentum. Thus the incoming b partons are given zero transverse momentum at leading order, and 
acquire transverse momentum at higher order. The two calculational schemes represent different pertur- 
bative expansions of the same physical process, and therefore should agree at sufficiently high order. It 
is satisfying that the NLO (and NNLO) calculations presented here agree within their uncertainties. This 
is a major advance over several years ago, when comparisons of bb — ► h at NLO and gg — > bbh at LO 
were hardly encouraging [1, 16]. 
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